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An impedance cardiograph (10) which determines cardiac output from a measurement of variations of chest impedance (2), provides 
method of calculating the effect of patient volume on the measurement from measurable patient height (H) and chest circumference (C), 
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IMPEDANCE CARDIOGRAPH APPARATUS AND METHOD 

Field of the Invention 
The present invention relates to impedance 
cardiographs which determine cardiac output by evaluating 
changes in impedance across the patient's chest cavity. 

5 Background of the Invention 

Impedance cardiography is a non- invasive technique 
for determining cardiac performance in humans. When such 
equipment is* employed, a high frequency electric signal 
is applied to the patient across outer electrodes 

10 positioned, for example, on the patient 1 s head and lower 
thorax. Voltage differences between sensing inner 
electrodes positioned between the outer electrodes on the 
patient's neck and chest are measured and used to compute 
an impedance (Z) . The impedance is based on the low 

15 magnitude, known electrical current passing between the 
outer electrodes. 

In 1932, Atzler and Ley-man reported that cardiac 
output of a human could be determined by such impedance 
methods in UJber ein neues Verfahren zur Darstellung der 

20 HerztatigJceit (Dielektrographie) , Arbeitsphysologie, 
5:636-680. In 1966, Kubichek reported the ability to 
correlate changes in base line impedance and the first 
derivative of impedance to stroke volume (SV) according 
to the following Equation (1) disclosed in U.S. Patent 

25 3,340,867 and in the publication: Development and 
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Evaluation of an Impedance Cardiac Output System, 
Aerospace Medicine, 37:1208-1212. 

0 dZ 
P* L2 *dt* L vet 

sv . ^ U, 

where : 

5 p is the resistivity of blood ; 

L is the spacing between the sensing electrodes; 

Zq is an average or baseline impedance; and 

a. 
dZ 

^£ is the magnitude of the peak negative value of 

the time derivative of the impedance Z for a period of 
10 time, typically a second. 

Cardiac output may be deduced from stroke volume by 
multiplying the latter times the heart rate. 

Although the Kubichek formula provides a value that 
correlated with cardiac output, the absolute accuracy of 
15 the method remained doubtful and, in particular, subjects 
with certain cardiovascular problems show values with 
great inaccuracies . 

In 1982, Sramek proposed a modification of the 
Kubichek formula of equation (1) which resolved the base 
20 line impedance Zq 2 into a dynamic and static component as 

reported in the publication: Cardiac Output by Electrical 
Impedance, Med. Elect., 2:274-290. The static term Zqs 

was described by the following equation: 

z 0 s - ^ «> 

25 where A is the area of the thorax being measured. 

The dynamic component ZQd was simply the baseline or 

average of the impedance being measured: 
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Zod = Z 0 (3) 

Incorporating the static term and dynamic term into 

the Kubichek equation provides the following formula: 

A. 

„ dZ 

sv = — r 0 — (4) 

5 The value of A may be estimated by approximating the 

chest as a cylinder in which case equation (4) becomes: 

A 

0 dZ 
c2 * L *dF* L vet 

SV = T^T Q (5) 

where C is circumference of the chest near the area 
of measurement. Alternatively, Sramek proposed that the 
10 term 

C 2 *L L 3 (1 7*H) 3 

-J^- be replaced with either J^S* or — ^25 — where H is 

the height of the patient because 1.7*H approximates L. 

These approximations did not produce good results 
and so Sramek was ultimately led to produce a set of 
15 charts attempting to establish correlation between area 
and the three factors of gender, height and weight. 

In 1986, Bernstein proposed a modified equation in 
which the separation of the electrodes and the height of 
the patient were considered, in the following form: 



[ l+1.5*h ]3 dZ 
I 2 J *dt* L vet 



20 sv = T^z~ 0 <s> 

All of the above methods suffer from lack of 
accuracy and indicate, in some subjects, falsely high or 
low values of stroke volume. 
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Summary of the Invention 
The present invention provides an improved method 
and apparatus for deducing stroke volume (and hence 
cardiac output) from impedance measurements. The 
5 invention provides an improved estimation of body volume 
and a processing of the derivative of the impedance 
signal that improves the reliability of the derived 
values of stroke volume and cardiac output . 

Specifically, in an impedance cardiograph for use on 

10 a human patient having a height of H and a chest 

circumference C, a means for applying an electrical 
excitation signal to the patient is used in conjunction 
with electrodes positioned on the chest with a separation 
distance of L, the electrodes producing a first 

15 electrical impedance signal Z which varies with impedance 
changes in the patient. A user input device, such as a 
keyboard, is provided to enter data on the height, 
electrode separation distance and circumference to 
produce corresponding second and third electrical signals 

20 H, L and C indicating those values. An electrical 

circuit, which may be an electronic computer, receives 
the second and third electrical signal and provides an 
indication of the patient's cardiac stroke volume, SV, as 
a function of Z, A, L and left ventricular ejection 

25 interval Lvet where A is deduced by the following 

approximation : 
L 

A " K*H*C 

where K is a predetermined constant. 
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Thus, the impedance cardiograph may calculate SV 
according to the following formula: 

L p*C 2 *L*L vet dZ' 
SV « K *H*C 4*7t*Z 0 dt 

which may be simplified to: 

_ P*C*L 2 *Lvet dZl 
5 SV " 4*K*H*7t*Z 0 dt 

or by combining constants: 

p»C*L2*Ly et d g, 
SV = K*H*Z 0 dt 

Thus, it is one object of the invention to provide a 
simple, yet more accurate characterization of a critical 

10 term A used in the calculation of stroke volume from 

patient impedance. In the present invention, the area of 
the impedance measurement, which is difficult to measure, 
is accurately derived from the readily measured values of 
chest circumference and patient height. 

15 The present invention has also recognized that 

variation in the derivative of impedance, one of the 

factors used in deducing stroke volume, is a significant 

source of inaccuracy in the computed stroke volume. 

dZ 1 

Accordingly, whereas may simply be the magnitude of 

20 the minimum derivative of impedance with time, it may 

also be compressed to reduce the amount that this maximum 
deviates from the norm. One method of weighting the 
maximum is according to the formula: 
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A. 

dZ 

dZ 1 dZn dt 



dt " dt 



where : 



/dZ.dZ 
dt'dt 



^ is the magnitude of the minimum time derivative 

Of Z; 

dZn . , dZ 
5 is the compressed value of 

dZ 

^ is a predetermined normal value and is an average 

A 

dZ 

value of ^£ for a population. 

Thus, it is yet another object of the invention to 
reduce the effect of a significant source of error in the 

10 calculation of stroke volume from impedance by 

implementation of a normal based weighting system. 

The foregoing and other objects and advantages of 
the invention will appear from the following description. 
In the description, reference is made to the accompanying 

15 drawings which form a part hereof and in which there is 
shown by way of illustration, a preferred embodiment of 
the invention. Such embodiment does not necessarily 
represent the full scope of the invention, however, and 
reference must be made therefore to the claims herein for 

20 interpreting the scope of the invention. 

Brief Description of the Drawings 
Fig. 1 is a block diagram of the circuitry for an 
impedance cardiograph according to the present invention, 
showing a computer as used to analyze the impedance data 
25 to produce a value of stroke volume and cardiac output; 
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and 

Fig. 2 is a flow chart of the software employed by 
the computer of Fig. 1 in analyzing the impedance and ECG 
signals acquired from the patient. 

5 Detailed Description of the Prefe rred Embodiment 

Referring to Fig. l, a portable impedance 
cardiograph 10 is connected to a patient 12 by five patch 
electrodes 14 (a) -(e). The first electrode 14(a) may be 
positioned on the patient's skin behind the right ear at 

10 the level of the ear canal. The second electrode 14(b) 
may be located at the left side of the neck on a flat 
surface approximately between the level of the chin and 
base of the hairline at least 5 cm below the level of 
electrode 14 (a) . The third electrode 14 (c) may be 

15 located just above the base of the sternum on the 

anterior median lines. The fourth and fifth electrodes 
14 (d) and 14 (e) may be located at least 5 cm below the 
electrode 14(c) on the patient's right and left sides 
respectively at the costal arch and the anterior axillary 

20 line. 

Electrodes 14 (a) and (d) are attached to an 
oscillator 18 which produces a constant current of 
approximately 1 milliamperes RMS through the patient 12. 
This electrical excitation establishes a series of equal 
25 potential surfaces through the patient 12 perpendicular 
to a line extending between the two outer electrodes 
14(a) and 14(d) . 

Electrodes 14(b) and 14(c) may sense the equal 
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potential lines generated by the current flowing between 
the outer electrodes 14 (a) and 14 (d) . Because the 
current between the outer electrodes 14(a) and 14(d) is 
of constant amplitude, the amplitude of the voltage 
5 sensed between the inner electrodes 14(b) and 14(c) is 
proportional to the thoracic impedance of the patient 12. 
The inner electrodes 14(b) and 14(c) are connected to 
differential amplifier 20 producing a signal Z. 

The amplifier 20 includes isolation circuitry that 

10 electrically isolates the inner electrodes 14 (b) and 
14(c) from the subsequent circuitry. Differential 
amplifier 20 also includes a precision half wave 
rectifier and low pass filter so as to provide a slowly 
varying DC signal whose value is proportional to the 

15 impedance being measured. The input impedance of the 

differential amplifier 20 is very high (e.g. 10 megohms) 
as compared to the impedance of the patient 12 between 
the inner electrodes 14(b) and 14(c). Thus, negligible 
current will flow through the inner electrodes 14 (b) and 

20 14(c) to amplifier 20. 

The impedance signal Z is received by a multiplexer 
32, such as are known in the art, to be periodically 
connected to an analog to digital converter 22 which 
samples the signal and provides a binary data word that 

25 may be read by microprocessor 24 via a bus 26. 

A first, vertical ECG signal is measured across 
electrodes 14(b) and 14(c) lying generally along a 
generally vertical line. This ECG signal is received by 
differential amplifier 30 to produce an 
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electrocardiograph signal E v according to techniques well 
known in the art . 

A second, horizontal ECG signal is measured across 
* electrodes 14(c) and 14(e) lying generally along a 

5 horizontal line. This ECG signal is received by 

♦ 

differential amplifier 31 to produce an 

electrocardiograph signal E^. Both signals E v and E^ are 
received by multiplexer 32 which periodically connects 
these signals to the analog to digital converter 22 to be 

10 sampled and converted to digital words for transmission 
on internal bus 26. 

Differential amplifiers 30 and 31 receiving the ECG 
signal also include isolation circuitry and a low pass 
filter having a cut-off frequency such as to 

15 substantially remove the 50 kHz oscillator signal from 
oscillator 18. As will be described further below, the 
direct and quadrature ECG signals are combined to produce 
a single ECG signal largely independent of the electrical 
orientation of the patient's heart, 

20 Also attached to bus 26 is computer memory 34 which 

may be composed of both random access memory ("RAM") and 
read only memory ("ROM") according to well known computer 
architectures. Memory 34 provides a means for storage of 
the binary representations of signals Z and E v and E^ 

25 under the control of microprocessor 24, and also holds a 
stored program defining the operation of the 
microprocessor 24 for the calculation of cardiac output, 
as will be described. 

Also attached to bus 26 is a pushbutton switch 25 
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which may be used by the patient to mark the occurrence 
of some event, such as a cardiac episode, during the 
recording of data from the patient 12 as will be 
described. 

5 The bus 26 also communicates with an infra-red 

transceiver 36, which permits the microprocessor to 
transmit and receive data to and from a similar 
transceiver 35 in a detachable base unit 37. The 
transceiver 35 is connected by a serial cable 38 to a 

10 desk-top computer 42 having a video monitor 43, a 

keyboard 44 and a disk drive 48 such as are known in the 
art. As such the impedance cardiograph 10 is portable 
and may be powered by internal batteries so as to be 
carried with the patient in the manner of a Holter 

15 monitor. 

Generally, the impedance cardiograph receives 
signals from the patient 12, isolates, amplifies and 
filters those signals, and then translates the signals to 
digital values which may be read and stored by 

20 microprocessor 24 to be processed according to the stored 
program in memory 34. Results of the processing may be 
transmitted to the computer 42 to be displayed on the 
video monitor 43 or saved on disk 48. 

The operation of the impedance cardiograph 10 

25 according to the stored program is controlled by a human 
operator through keyboard 44. The operator prepares the 
patient 12 for the impedance measurements and may enter 

■ 

certain data to the keyboard 44 that characterizes the 
patient 12 and that is necessary for the analysis of the 
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signals from the patients 12 as will be described. This 
analyses is done in partially in the microprocessor 24 so 
as to reduce the amount of data to be stored in memory 
34, but, as will be understood in the art, the analyses 
5 of the data may be shared between the microprocessor 24 
and the computer 42 as a matter of engineering choice. 

Referring now to Fig. 2, at the first step in the 
analyses, indicated by process block 50, certain data 
related to the particular patient 12 or related to 

10 fundamental and essentially universal physiological 
parameters, may be entered by the operator. 

The first of these parameters is p which is the 
resistivity of blood in ohms-cm. Generally, this value 
may be approximated as a constant for all patients, 

15 however, it may be modified by the operator in abnormal 
cases based on the measurement of hematacrit . 

A second value, C, is the circumference of the 
patient f s chest in cm taken at the site of inner 
electrode 14(c) around the patient's chest. 

20 The value L is also input, being the distance 

between the inner electrodes 14(b) and 14(c) in cm. The 
height of the patient, H, in cm is also entered. 

Once the necessary fixed parameters are entered at 
process block 50, two variable parameters: heart rate HR 

25 and ventricular ejection time L ve t are entered per 

process block 52. For this purpose, the ECG signal E may 
be directly displayed on the video monitor 43 so that 
these quantities may be determined according to methods 
well known in the art. The ECG signal is calculated from 
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the vector sum of the values of E v and E^ most simply as 
follows: 

ECG = + Eg 

This vector summing reduces the need to precisely 
5 orient the ECG electrodes with respect to an electrical ♦ 
polarity of the heart and therefore in practice provides 
a superior ECG signal. 

Generally, Lvet ^ s t * ie time between the opening of 

the aortic valve and the closing of the aortic valve. 

10 The heart rate is simply the number of beats per second 
which is the inverse of the period between successive R 
waves . The heart rate may be averaged over a number of 
beats according to methods well known in the art. Both 
quantities may be determined by inspection by the 

15 operator or preferably may be determined automatically 
after sufficient ECG and Z data is acquired as indicated 
by process block 52 . The value of Lvet in the preferred 

embodiment is determined by analyzing the impedance 

dZ 

signal Z to measure a period beginning when ^£ is first 

dZ 

20 less than zero and ending when the value of ^ reaches a 

local maximum above zero. The heart rate HR is measured 
by detecting and counting R waves in the ECG signal. 

The acquisition of the ECG and Z signal per process 
block 54 continues. This acquisition is on a continuous 
25 basis and occurs concurrently with the subsequent 
calculations so that the cardiac output may be 

c 

continuously stored in memory 34 or displayed in 
essentially a real-time manner. 
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The acquired impedance data is in the form of 

discreet samples taken approximately 300 times per 
second, each sample which may be represented by Z± where 

'» i is an index number of the particular sample. As each 

5 sample is acquired, it is stored in consecutive 

addresses in memory 34 to indicate its relative position 
with respect to other samples and to indicate the time of 
the sample indirectly through the constant sampling rate. 
Because the impedance cardiograph 10 is portable, 

10 there is a risk that artifacts may be introduced into the 
impedance measurement by electrode movement. This is 
because the measured impedance values are approximately 
two orders of magnitude lower than the electrode to skin 
resistance. Accordingly, the impedance data over a 

15 period of approximately one minute is "ensemble" 

averaged, per process block 56, to reduce its noise 
content. Ensemble averaging is a well known technique in 
which blocks of impedance data are averaged on a point - 
by-point basis with other blocks of impedance data so 

20 that the averaged points are from corresponding portions 
of the impedance waveform cycle. Thus, the shape of the 
impedance waveform is not destroyed in the averaging 
process. In order to perform such ensemble averaging, it 
is necessary to identify a common fiducial point to align 

25 the blocks of data. Selection of the fiducial point must 
be extremely precise, otherwise the characteristics of 
the impedance waveform will be "blurred" by a mis- 
registration of other blocks. 

This fiducial point may be the peak of the R wave of 
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the ECG signal. Normal techniques for determining the 
time of the R wave, such as may be used for the 
measurement of heart rate, however, are not suitably 
accurate for the purpose of ensemble averaging. r 
5 Accordingly an extremely accurate identification process 
is used. First, as represented by process block 58, the 
ECG signal is monitored to isolate a standard R wave. 
Only portions of the received ECG signal having no 
detectable artifacts or noise are considered. This 

10 standard R wave is then correlated to the incoming ECG 

signal to identify the precise location of the R wave (by 
the value of highest correlation) . This location is used 
at the point of common alignment for the impedance 
waveforms to be ensemble averaged. Periodically, a new 

15 standard R wave is obtained so that the standard remains 

current over time. 

After pairs of data Zi and Zi+1 are acquired and 

dZi 

averaged, a derivative value is may be computed by a 

simple subtraction of adjacent samples of the ensemble 
20 average per process block 58, that is: 
dZi 

"dT = Z i+1 " Z i (7) 
Alternatively, in order to reduce the presence of 50 

Hz or 60 Hz noise, this derivative computation can employ 

samples Zi+g - Zi or Z±+5 f respectively. 

25 It has been determined that variation in the 

A 

dZ . 

magnitude of the minimum of this value, is a 

significant source of error in the calculation of cardiac 
output. Accordingly, at process block 60, a compressed 
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dZn 

derivative, is computed according to the following 

formula : 



dZ 

dZn dt 
dt = 



(8) 



fdZ dZ 
dt'dt 



where : 



dZ 

5 dt ^ s t * le magnitude of the minimum time derivative 

of Z as previously defined; 

dZn dZ 
is the compressed value of 

dZ 

^ is the predetermined normal maximum value and is 
dZ 

an average value for a population and is about 1.73 

10 ohms per second. 

Equation (8) has the effect of reducing the 
dZ dZ 

excursions of ^ from Other nonlinear compression 

systems may also be used provided they have the effect of 
dZ ^ 

compressing ^£ about the norm. 

15 At succeeding process block 52, the value of Zq is 

also computed. Zq is the base transthoratic impedance 
and in this implementation, simply the average value of Z 
for one cardiac cycle . Because of the need to average a 
number of samples, when the sampling is first begun, no 

20 display is provided to the video terminal 43 until 

sufficient samples have been made to insure the accuracy 
of this value Zq. 

At process block 64 a stroke volume may be 
calculated according to the following formula: 
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L P*c2*L*Ly et dZ 
SV ~ K*H*C 4*7t*Z 0 dt {9) 

It will be recognized that this is simply Equation 
(5) with the addition of a factor p and the addition of a 
factor k*h*C' t ^ 1 ^ s l atter factor approximating a slice of 

5 body volume in the area of the impedance measurement. 
This calculation of stroke volume has shown 
significant improvements in accuracy and correlation 
coefficients in clinical studies. 

Also at process block 64, cardiac output may be 
10 determined by multiplying the stroke volume times the 
heart rate: 

CO = SV*HR (10) 
As cardiac output is computed, it is displayed in 
graphical form on video monitor 43. Thus, the operator 
15 is provided with concurrent ECG data and cardiac output 
data on an essentially real-time basis per the display 
indicated by process block 66. 

After each updating of the display of 43 is 
accomplished, the program acquire additional data until 
20 the measurement session is complete. 

While this invention has been described with 
reference to particular embodiments and examples, other 
modifications and variations will occur to those skilled 
in the art in view of the above teachings. Accordingly, 
25 the present invention is not limited to the preferred 
embodiment described herein, but is instead defined in 
the following claims. 
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Claims 

I claim: 

1. An impedance cardiograph used to evaluate 
cardiac output of a human patient having a height and a 
chest circumference and chest cross-sectional area A 
comprising: 

means for applying an electrical excitation signal 
to the chest of a patient; 

electrodes, adapted to be positioned on the patient 
with a separation distance of L thereby being responsive 
to the excitation signal, for producing a first 
electrical signal Z which varies with impedance changes 
in the patient; 

input device for receiving values of the height and 
circumference and providing corresponding second and 
third electrical signals H and C indicating those values; 

electrical circuit means for receiving the first, 
second and third electrical signal and providing an 
indication of the patient's cardiac stroke volume SV as a 
function of Z, A, and L; 

where A is deduced by the following approximation: 
L 

K*H*C 

where K is a predetermined constant. 
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2 . The impedance cardiograph of claim 1 wherein the 
electric circuit means calculates SV according to the 
following formula: 

L P*c2*L*Lv et d g, 



SV = 



K*H*C 4*n*Z n dt 



where 

p is the resistivity of blood; 

Zq is an average baseline impedance; 

Lvet is a left ventricular ejection interval; and 
dZ' 

is a function of the time derivative of the 



impedance Z. 



dZ' 

3 . The impedance cardiograph of claim 2 wherein 
is the magnitude of the minimum time derivative of Z. 



dZ' 

4 . The impedance cardiograph of claim 2 wherein 
is a compressed value of a time derivative of Z about a 

A. 

dZ 1 

predetermined normal value so that has a value less 
than the time derivative of Z if the time derivative of Z 

A 

dZ' 

is greater than a predetermined normal value and has 

a value greater than the time derivative of Z if the time 
derivative of Z is less than the predetermined normal 
value . 
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5 . The impedance cardiograph of claim 1 wherein K 
is substantially 0.875. 



6 . An impedance cardiograph used to evaluate 
cardiac output of a human patient comprising: 

means for applying an electrical excitation signal 
to the chest of a patient; 
5 electrodes, adapted to be positioned on the patient 

for producing a first electrical signal 2 which varies 
with impedance changes in the patient; 

a non-linear processor for calculating a time 
derivative of Z, and providing a compressed derivative, 
10 having a value less than the time derivative of Z if the 
time derivative of Z is greater than a predetermined 
normal value and having a value greater than the time 
derivative of Z if the time derivative of Z is less than 
the predetermined normal value; 
15 electrical circuit means for providing an indication 

of patient's cardiac stroke volume SV as a function of 
the compressed time derivative. 
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7. The impedance cardiograph of claim 6 wherein the 

non- linear processor produces the compressed value of the 

time derivative of Z according to the following formula: 

dZ 

dZn dt 



dZ dZ 
dt'dt 



dt = 
where : 

A 

dZ 

^ is a time derivative of Z 
dZn 

is the compressed time derivative of Z 

dZ 

is the predetermined normal value and is an 

A. 

dZ 

average value of for a population. 



8 . The impedance cardiograph of claim 6 wherein the 
predetermined normal value is substantially 1.73 ohms per 
second. 
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9. A method of evaluating the cardiac output of a 
human patient having a height H and a chest circumference 
C and chest cross -sectional area A comprising the step 
of: 

(a) applying an electrical excitation signal to the 
chest of a patient; 

(b) positioning electrodes on the patient with a 
separation distance of L thereby being responsive to the 
excitation signal to produce a first electrical signal Z 
which varies with impedance changes in the patient; 

(c) employing an electronic computer to: 

(i) receive values of the height and 
circumference and providing corresponding second and 
third electrical signals H and C indicating those values; 

(ii) evaluate the patient's cardiac stroke 
volume SV as a function of Z, A, L and left ventricular 
ejection interval L V et; 

where A is deduced by the following approximation: 
L 

K*H*C 

in which K is a predetermined constant. 



WO 96/01586 PCT/US95/08856 

-22- 

10 . The method claim 9 wherein the electronic 
computer calculates SV according to the following 
formula : 

L P*c2*L*L vet dZ' 
SV = K*H*C 4*7t*Z 0 dt 

where 

p is the resistivity of blood; 

Zq is an average baseline impedance; and 

dZ* 

is a function of the time derivative of the 
impedance Z. 

dZ' 

11. The method of claim 10 wherein is the 
magnitude of the minimum time derivative of Z. 



12. The method of claim 10 wherein is a 

compressed value of the time derivative of Z about a 

dZ' 

predetermined normal value so that has a value less 

than the time derivative of Z if the time derivative of Z 

dZ' 

5 is greater than a predetermined normal value and ^r - has 

a value greater than the time derivative of Z if the time 
derivative of Z is leiss than the predetermined normal 
value. 
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13 . The method of claim 12 wherein K is 
substantially 0.875. 



14. A method of evaluating cardiac output of a human 
patient comprising the steps of: 

(a) applying an electrical excitation signal to the 
chest of a patient; 

(b) positioning electrodes on the patient for 
producing a first electrical signal Z which varies with 
impedance changes in the patient; 

(c) determining a compressed derivative of the 
impedance signal Z, having a value less than the time 
derivative of Z if the time derivative of Z is greater 
than a predetermined normal value and having a value 
greater than the time derivative of Z if the time 
derivative of Z is less than the predetermined normal 
value ; and 

(e) evaluating the patient's cardiac stroke volume 
SV as a function of the compressed time derivative. 
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15. The method of claim 14 step (c) is according to 
the following formula: 

-A. 

dz 

dZn dt 



dt -\ /d| ,dZ 
V dt' dt 



5 where : 

A 

dZ 

^ is the magnitude of the minimum time derivative 

of Z 

dZn , 

is the compressed time derivative of z 

dz 

^ is the predetermined normal value and is an 

.a 
dZ 

10 average value of ^£ for a population. 



16. The method of claim 14 wherein the 
predetermined normal value is substantially 1.73 ohms per 
second. 
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